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ABSTRACT 

We investigate the dispersion in the internal colours of faint spheroidals in the 
Hubble Deep Fields North and South. In high redshift rapid-collapse scenarios, 
the dispersion in internal colours should be small at moderate redshift apart from 
a small metallicity-induced reddening in the enriched cores. However, recently- 
assembled spheroidals are likely to show non-homologous internal colours, at least 
until younger stellar populations become fully mixed. Here we find that a re- 
markably large fraction ( ^ 30%) of the morphologically-classified spheroidals with 
huw < 24 mag show strong variations in internal colour, which we take as evidence 
for recent episodes of star-formation. In most cases these colour variations manifest 
themselves via the presence of blue cores, an effect of opposite sign to that expected 
from metallicity gradients. Examining similarly-selected ellipticals in five rich clus- 
ters with 0.37 < z < 0.83 we find a significant lower dispersion in their internal 
colours. This suggests that the colour inhomogeneities have a strong environmental 
dependence being weakest in dense environments where spheroidal formation was 
presumably accelerated at early times. We use the trends defined by the cluster sam- 
ple to define an empirical model based on a high-redshift of formation and estimate 
that at z ~ 1 about half the field spheroidals must be undergoing recent episodes of 
star-formation. Using spectral synthesis models, we construct the time dependence 
of the density of star-formation (psfr)- Although the samples are currently small, we 
find evidence for an increase in psfr between z = to z = 1 . We discuss the implica- 
tions of this rise in the context of that observed in the similar rise in the abundance of 
galaxies with irregular morphology. Regardless of whether there is a connection our 
results provide strong evidence for the continued formation of field spheroidals over 
0< z <1. 
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1 INTRODUCTION 

The age distribution of elliptical galaxies remains a central 
topic in galaxy evolution. For many years ellipticals were 
viewed as old systems which formed as the result of a short 
and intense burst of star formation (Baade 1957, Sandage 
1986), after which the stars within them passively evolved. 
Strong observational support for this view came from old 
and coeval ellipticals found in rich clusters at low (Bower 
Lucey & Ellis 1992) and high redshifts (Ellis et al. 1997, Stan- 
ford et al. 1997). 



* fmenante@ast.cam.ac.uk 



However, in hierarchical models of galaxy formation 
dominated by cold dark matter (CDM), elliptical galaxies 
form over a longer period, as the result of mergers between 
low mass disks (Kauffmann et al. 1996, Baugh et al. 1996). 
In these models clusters form from the highest peaks of the 
density fluctuations, so the existence of old, coeval cluster el- 
lipticals is not necessarily in contradiction with hierarchical 
models. Clearly, regions of high density are poorly suited for 
testing models for elliptical galaxy formation. 

The controversy regarding the epoch of formation of 
elliptical galaxies has prompted several observational cam- 
paigns designed to establish whether field ellipticals share 
the same evolutionary history as their clustered counter- 
parts. Recent tests have focused on the very red optical- 
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IR colours predicted by single-collapse models. Zepf (1997) 
and Barger et al. (1998) adopted a colour-based approach 
by searching for a very red tail in the optical-IR colour 
distribution of faint HDF galaxies. Very few sources were 
found with colours matching those expected for high red- 
shifts of formation. Similarly Menanteau et al. (1999), using 
a statistically-complete sample of morphologically-selected 
spheroidals from Hubble Space Telescope (HST) archival 
data, concluded on the basis of optical-infrared colours that 
field ellipticals cannot have formed the bulk of their stars in 
a single-burst of formation at a high redshift. 

However, an early period of collapse may be rescued 
if massive systems subsequently suffer minor episodes of 
star formation rendering the observed colours bluer than the 
limits explored in the above observational studies. Jimenez 
et al. (1998) propose a multi-zone model of spheroidal for- 
mation which predicts bluer colours over the redshift range 
explored. This is consistent with recent studies by Abra- 
ham et al. (1999) and Kodama, Bower & Bell (1998, hereafter 
KBB) which found a large fraction (~ 30%) of spheroidals 
have properties which differ from those predicted by single- 
collapse high-redshift models. More recently Tamura et al. 
(2000) studied the origin of colour gradients in bright ellipti- 
cal galaxies in the Northern HDF where they report a simi- 
lar fraction (~ 30%) of systems with blue integrated colours. 
However, they claim that colour gradients are very small and 
originated by stellar metallicity. 

The ability of HST to resolve distant galaxies, and the 
very deep multi-colour photometry available in the Hubble 
Deep Fields (HDFs), opens up new possibilities for probing 
the internal characteristics of distant galaxies. In this paper 
we adopt a methodology similar to that used in Abraham 
et al. (1999), where the spatially-resolved colours of galax- 
ies of known redshift in the Northern HDF were used to 
probe their star-formation history. We improve on Abraham 
et al. (1999) by comparing the internal properties of field 
spheroidal with their clustered counterparts and by enlarg- 
ing the sample of early type galaxies by a factor of 6, includ- 
ing the incorporation of high signal-to-noise (S/N) data from 
the Southern HDF. 

A plan for the paper follows. In §^ we present the field 
and cluster spheroidal sample. In §^ we introduce the prin- 
ciples of our methodology and discuss our results and their 
implications for the evolution of spheroidals under a free- 
model approach. In §[j| we attempt to model the observed 
colour dispersions, estimate the corresponding time-scales 
for the star-formation activity and compute the present SFR 
for spheroidals from the HDFs. In §^ we summarize our con- 
clusions. 



2 FIELD AND CLUSTER SAMPLES 

2.1 Spheroidals in the Hubble Deep Fields 

The Hubble Deep Fields (HDFs) provide our primary sam- 
ple of field faint spheroidals. In the present paper sources in 
the HDF fields are designated using the IAU IDs given in 
version 2 of the HDF SExtractor catalogues constructed 
by the Space Telescope Science Institute. Vega magnitudes 
are used throughout this paper. HDF galaxies were morpho- 
logically classified using both the automated methodology 



Table 1. The HST Cluster Sample 



Cluster 


Filter 1 


Filter 2 


z 


A370 


F814W 


F555W 


0.37 


C10412-65 


F814W 


F555W 


0.51 


C10016+16 


F814W 


F555W 


0.55 


C10054-27 


F814W 


F555W 


0.56 


MS1054+03 


F814W 


F606W 


0.83 



described in Menanteau et al. (1999) (based on both, the cen- 
tral concentration (C) and asymmetry (A) parameters from 
Abraham et al. 1996), as well as using the visual classifica- 
tions made by one of us (RSE). Visual and automated clas- 
sification have been shown to agree well in previous studies 
(eg. Brinchmann et al. 1998; Menanteau et al. 1999). How- 
ever, Marleau and Simard (1998) claim that a significant frac- 
tion of galaxies visually classified as ellipticals in the van der 
Bergh et al. (1996) morphological catalog of the HDF were 
disk-dominated galaxies, although most of these systems are 
fainter than our selection limits. 

A limiting magnitude of I ni = 24 mag was adopted, so 
that information for each galaxy is distributed over a large 
enough number of pixels to allow us to define meaningful 
measures of our model-independent estimator that we will 
describe in section ^ We augment the spectroscopic red- 
shifts in our sample using photometric redshift estimates 
kindly provided by S. Gwyn (Gwyn et al. 2000 in prepa- 
ration). Photometric redshifts play a particularly important 
role in the Southern HDF field, where the published spec- 
troscopic data is presently rather limited. Our final sample 
consists of 79 spheroidal galaxies, 24 of which have spectro- 
scopic redshifts and 55 of which have photometric redshifts. 

2.2 Cluster Spheroidals 

In order to compare field and cluster samples directly, we se- 
lected a sample of high-redshift cluster ellipticals observed 
in two passbands with HST's Wide Field Planetary Camera 
2 (WFPC2) over a similar wavelength baseline to V 606 — J sl4 
and spanning a redshift range comparable to the one ob- 
served for fi eld spheroidals. F rom the MORPHS collabora- 
tion dataset (Smail et al. 1997), we selected 4 rich cluster at 
0.37 < z < 0.56. We also used the data from van Dokkum 
(1999) for the X-ray cluster MS1054+03 at z = 0.83. 

Both studies have published morphological classifica- 
tions, from which we select spheroidals as objects classified 
as E and E/S0 (but not SO). In Table [3] we list the properties 
of the selected clusters. It is important to note that the cluster 
imaging data is not as deep as the HDF data, and to allow for 
this our cluster sample is restricted to be brighter than a limit- 
ing magnitude of J S i4 = 22 mag, ie. two magnitudes brighter 
than the corresponding HDF field data. We discuss the po- 
tential limitations introduced by adopting different cluster 
and field limiting magnitudes further below. Another differ- 
ence is that the MORPHS sample was imaged using F555W 
instead of the F606W filter which defines the V-band in the 
HDF. We defer a detailed study of its influence on our analy- 
sis until the next section. 
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Figure 1. The methodology of internal colour variations as applied to two 
HDF-N spheroidals. The upper panel shows an Igu = 20.48 mag example 
with low internal scatter at z — 0.77, whereas the bottom panel illustrates 
an Igi4 = 21.66 example with a bluer core at z = 0.48. The pixel-by-pixel 
colour distributions are shown alongside each case. Coloured dots refer to 
pixels whose SNR is above the adopted threshold. 



3 MODEL-INDEPENDENT ANALYSIS 
3.1 Methodology 

In some ways our methodology for examining the internal 
colours of an individual galaxy can be considered as a gener- 
alization of the use of colour-magnitude relations in studying 
the history of ellipticals in galaxy clusters: the photometric 
dispersion is a powerful probe of variations in star-formation 
history provided systematic effects are under control (Bower, 
Lucey & Ellis 1992). 

Initially we will focus on a largely model-independent 
approach to the problem of understanding how the inter- 
nal colours of field ellipticals vary with redshift relative to 
those of their clustered counterparts. We will concentrate 
only on J sl4 and V 606 bands, as the resulting V 606 — h\n 
colours have substantially smaller observational errors than 
the B 450 — V 606 and L7 30 o — B450 colours, especially for systems 
dominated by old stellar populations. 



3.1.1 The Homogeneity Estimator 

After registering the Northern HDF images to sub-pixel ac- 
curacy we isolated each galaxy from the background sky by 
selecting all contiguous pixels within a surface brightness 
threshold of V 606 = 25.5 mag arcsec~ 2 . We select this limit in 
order to maximize the signal/noise per pixel associated to 
V 60 6 - hu colours. 

The distribution of V 60 e — hu colours for individual 
pixels in the galaxy is used to calculate a quantity 5(V — I) 








Figure 2. The observed S(V — I) as a function of redshift for modelled 
galaxies viewed at various redshifts with two photometric colour systems. 
The V606 — Jsi4 (solid lines) and V555 — Ign (dashed line) curves are plot- 
ted for two cases: a single-burst ("Passive") system corresponding to an 
epoch of formation Zp = 5, and an "Active" case which incorporates a 15% 
by mass burst of star formation centrally super-imposed on an old galaxy 
at the time of observation. 



which characterises the internal homogeneity of a galaxy. A 
5(V — I) statistic is defined as: 

Z(x t - x) 2 S(xi)SNR{xi) 



5(V-I) =2N 



ZS{xi)SNR{xi) 



ZxiS(xi)SNR(xi) 
Y.S(xi)SNR{xi) 



(1) 



(2) 



where S(x) is a selection function for pixels with sig- 
nal/noise above a certain threshold (SNR > 1.3) such that 
S(x) = 1 for pixels above and S(x) = for pixel below 
the threshold. SNR(xi) is the signal-to-noise ratio for a given 
pixel colour x u and N is an arbitrary scale factor. The selec- 
tion function S(x) and the weighting according to SNR(x) 
address biases arising from noise variations at pixel scales by 
rejecting low signal pixels and weighting the pixels contribu- 
tion proportionally to their signal. Figure [l] show an example 
of spheroidals in HDF-N with both high and low S(V — I) . 

As mentioned previously, the MORPHS dataset was ob- 
served with a different filter baseline than that used for the 
HDF. To quantify the effect of this difference on S(V — I), we 
used stellar population synthesis libraries (Bruzual & Char- 
lot 1996, BC96) to calculate the observed colours (in both 
photometric systems) for two galaxy types: a) an old single- 
burst elliptical (formation epoch z f = 5 and e-folding time 
t = 1 Gyr) and, b) a star-forming elliptical galaxy consisting 
of a recent burst (100 Myr) centrally superimposed on an old 
system. In order to emulate the geometrical properties of the 
modelled galaxies, we use the observed colours at a given 
redshift to create elliptical galaxies with r 1 / 4 profiles using 
a customised version of the IRAF package ARTDATA and a 
fixed physical size with a half-light radius r e = 2.5 kpc. For 
the case of the blue core, we attempt to match the properties 
of Figure [i] lower panel, with a burst that accounts for 15% of 
the total mass of the galaxy with and a physical size of 0.1r e . 
5(V — I) for both photometric systems are shown as a func- 
tion of redshift in Figure ^. This exercise demonstrates that 
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Figure 3. Tfe colour inhomogeneity 5(V — I) for real and simulated data 
sets as a function of their central concentration indicating the level of bias 
introduced by the wavelength-dependence of the HST point spread function. 



the substitution of F555W for F606W does not have a major 
impact on the calculation of 5(V — I) . The slightly longer 
wavelength baseline in the cluster 5(V — I) estimates should 
result in a greater sensitivity to recent activity. 

It is important to note that because 5(V — I) is based on 
observed (rather than rest-frame) colours, its sensitivity to star 
formation is redshift-dependent. Unfortunately, it is not pos- 
sible to use additional filters to create a rest-frame equiva- 
lent since the shorter wavelength HDF data has poorer sig- 
nal/noise and the cluster data only consists of two pass- 
bands. In order to quantify this restriction, we studied val- 
ues of 5(V — I) obtained for the blue core case modelled in 
Figure ^[ We placed this modelled galaxy, at the epoch corre- 
sponding to its maximum inhomogeneity (representing the 
most extreme case we will discuss) at various redshifts cal- 
culating 5(V — I) at each redshift. Figure ^ shows increased 
sensitivity to the fixed star formation rate at z ~ 0.5 for 
F606W and z ~ 0.4 for F555W; otherwise the sensitivity is 
broadly constant. This increase arises as the 4000A break 
in the underlying passive system passes through the F606W 
and F555W filters and will not seriously compromise our 
analysis. 



3.1.2 The Influence of the HST Point-Spread Function 

As the HST point-spread function (PSF) varies as a func- 
tion of wavelength, spurious centrally-concentrated inhomo- 
geneities might be expected for sharply-peaked profiles such 
as those encountered in spheroidal galaxies. In order to test 
this, we performed extensive Monte Carlo simulations using 
the IRAF package ARTDATA. Artificial galaxy images were 
created and subsequently analysed using procedures iden- 
tical to those used to analyse our observed data. 

Artificial spheroidal galaxies were synthesised with de 
Vaucouleurs profiles, using stars from the HDF-N as PSF 
templates. Simulated images were generated in both I$u and 
V 606 -bands. We explored a range of half-light radius (r e ) 
from 0.04" - 0.65" corresponding to ~ 0.4-6 kpc (Ho = 
65 km s _1 Mpc~ 3 , qa = 0.1), and a range of magnitudes rep- 



resentative of our sample (18 < J 8 i4 < 24 mag). As the in- 
fluence of the PSF on the central colours is expected to be 
a strong function of the "peakiness" of the central portion of 
the galaxy, we chose to probe the influence of the PSF on both 
real and simulated galaxies as a function of the central con- 
centration (C) parameter defined in Abraham et al. (1994). 

Figure ^ shows the result of this exercise. The PSF varia- 
tions affects the 5(V — I) statistic only for the most centrally 
concentrated (essentially stellar) objects. Over the observed 
range of C for the real data, 0.3 < C < 0.7, 5(V - I) in the 
simulated data is well below the values recovered from the 
real observations, and we conclude that this effect will not 
seriously affect any of the following conclusions. 

3.2 Evolution in Field Spheroidals 

Figure ^ shows V 606 — J 814 colour maps for the 79 spheroidals 
in our HDF sample sorted in ascending order according to 
Igi4 magnitude. Redshift and 5(V — I) values are also marked 
in each case; redshifts in parenthesis refer to photometric 
estimates. The most striking characteristic of these data is 
the large proportion with significant internal colour varia- 
tions. We emphasise that we have not significantly varied 
the "colour stretch" from panel to panel (as indicated in the 
colour bar) to exaggerate this result. Adopting, for simplic- 
ity, a colour variation of ~ 0.2 mag as indicative of the effect 
(since this broadly corresponds to a change in visual colour 
according to our look-up table), it appears that ~ 40% of the 
field spheroidals show a detectable degree of inhomogeneity. 
In most cases where inhomogeneities are seen, these appear 
to be in the cores. Obviously, the visual inhomogeneities cor- 
relate closely with 5(V — I) . 

In § 3.3 we will calibrate these inhomogeneities with ref- 
erence to a corresponding sample of cluster spheroidals, but 
it is meanwhile interesting to consider whether these inter- 
nal colour variations correlate with independent constraints. 
In doing this we will make the assumption that the bluer 
colours seen are indicative of recent star formation. 



3.2.1 Comparison with Kodama et al. 

We first compare our 5(V — I) estimates with the results 
of Kodama, Bower & and Bell (1998, KBB). Those authors 
used stellar population models and photometric redshifts 
to construct a rest-frame U — V vs. V colour-magnitude di- 
agram for field spheroidals in the HDF-N. They took the 
morphologically-classified sample of 35 spheroidals with 
K < 20.5 from Franceschini et al. (1998) of which 29 are in 
our sample. The 6 missing cases correspond to two classes 
of objects: a) objects that only Franceschini et al. classified 
as spheroidals, and b) double detections in the Frances- 
chini et al. (1998) study which we have eliminated. Using 
spectral synthesis models, KBB concluded that the sample 
could be described as comprising a red sequence of old 
passively-evolving systems, with an additional component 
(comprising 30% of the total sample) with bluer optical- 
infrared colours that lay significantly off the fiducial colour- 
magnitude relation. 

We have divided the Franceschini et al. subsample of our 
HDF-N spheroidals into 'red' and 'blue' examples depend- 
ing on whether their rest-frame colours lie within the red se- 
quence limits defined by KBB. Specifically, we classed as blue 
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Figure 4. Colour images of 79 HDF field spheroidals keyed to their Vgo6 — hu colours. The integrated 7gi4 mag for each galaxy is shown in each sub- 
together with the redshift and 5(V — I) . Redshifts in parenthesis indicate photometric estimates. Axes labels correspond to arcseconds 
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available at http://www.ast.cam.ac.uk/~fmenante/HDFs/mos_all_B.ps.gz 



Figure 4 - continued 



all spheroidals whose rest-frame colours have U — V < 0.8. 
Figure U compares our 5(V — I) with the red /blue classifica- 
tions demonstrating a strong correlation with the estimated 
U — V rest-frame colour from KBB. In all but one case, red 
KBB spheroidals correspond to low 5(V — I) systems in our 
scheme. We estimate that ~ 9 out of of 29, or ~ 30% of 
spheroidals have values of 5(V — I) suggesting some degree 
of star-formation activity, a value in excellent agreement with 
the KBB estimate. 

The strong correlation between our model-independent 
measures of 6(V — I) and KBB's estimates of rest-frame 
U — V is reassuring evidence that simple spectral synthesis 
models used to describe spheroidal evolution offers a valid 



way to quantitatively interpret our results. Armed with this 
reassurance, in §4.2 we will use similar models to attempt to 
constrain the amount of star-formation contributed by evolv- 
ing spheroidals as a function of redshift. 



3.2.2 Low 6(V — I) field spheroidals 

Despite the significant fraction of spheroidals showing some 
degree of colour variation, the majority show homogeneous 
internal colours, suggesting coeval evolution of their stel- 
lar populations. We selected the galaxies with low 5(V — I) 
< 0.04 and compared their integrated colours against the pre- 
dictions of a monolithic collapse model. As expected this 
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Figure 5. The distribution of 5(V — I) with redshift for those HDF 
spheroidals studied by Kodama et al. (1998, KBB). The size of the symbols 
indicates the l%uw integrated magnitude. 
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Figure 7. Mean &(V — I) values and Iff dispersions calculated for 
spheroidal galaxies in 5 clusters observed with HST in similar passbands to 
those used for the HDF field study. 
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Figure 6. The V$o6 ~ hu observed colours for a subsample of internally 
homogeneous (5(V — I) < 0.04) spheroidals as a function of redshift com- 
pared with the predictions of a high-redshift single burst model (zp = 5, 
T = 1 Gyr) and solar metallicity for a set of cosmologies. Inverted trian- 
gles refer to HDF-S and unfilled symbols refer to galaxies for which only 
photometric redshifts are available. 



sample closely tracks the expected colour-redshift relation of 
a passively evolving (Figure |^) providing reassurance that 
5(V — I) is a good statistic from which to differentiate active 
and passive systems. 

An interesting feature is the somewhat larger fraction of 
homogeneous spheroidals within HDF-N with respect to the 
HDF-S. Only 4 spheroidals in the HDF-S show low 5(V - I) 
values compared to 15 spheroidals in the HDF-N. This could 
be interpreted as evidence for clustering in the HDF-N as 
suggested by Kodama et al. (1998). 



3.3 Evolution in Cluster Spheroidals 

As we discussed in §1, it is now commonly-believed that 
spheroidals in the cores of rich clusters formed the bulk of 



their stars before z > 2 (Bower et al. 1992, Ellis et al. 1997, 
Governato et al. 1998, Kauffmann et al. 1999). Cluster sam- 
ples therefore represent suitable benchmarks whose 6(V — 1) 
values at a given redshift can be used to provide a model- 
independent calibration of the internal colour scatter ex- 
pected from an old stellar population at a given redshift. 

As mentioned earlier, our limiting magnitude in the 
cluster sample was chosen as Igu = 22, considerably brighter 
than adopted for the deeper HDF images. The morphological 
classifications were taken from the published articles. In the 
case of the MORPHS sample there is a cross-check given one 
of us (RSE) classified both samples. As in the previous section 
we have calculated 5(V — I) for each cluster spheroidal using 
the prescription given in equation ^[ We find a striking ab- 
sence of internally inhomogeneous systems. In fact, we find 
no cluster spheroidals with blue cores in any of the five clus- 
ters. We then calculate a median V 606 — 7 814 value and will 
take this as representative of the 5(V — I) value appropri- 
ate for an old spheroidal population over the redshift range 
0.37 < z < 0.83 (Figure |). 

The strongest result of this paper is that field spheroidals 
show substantially greater internal colour dispersions than 
their clustered counterparts. Since the limiting magnitude 
of the field sample is two magnitudes brighter than that 
adopted for the cluster sample, this is most cleanly demon- 
strated by comparing 5(V — I) values on a galaxy -by-galaxy 
basis in both samples as a function of rest-frame luminos- 
ity (Figure |). Absolute magnitudes were calculated assum- 
ing H = 65 km s _1 Mpc -1 . In the cluster analysis, galaxies 
smaller than 60 pixels were not considered, as the 5(V — I) 
measurement was determined to be too noisy to be mean- 
ingful. Although about 15% of the cluster spheroidals were 
rejected on this basis, from Figure ^ we are reassured that we 
are probing a similar range of Mg45o absolute magnitudes. 

For the cluster spheroidals we calculated Mg45o for each 
galaxy from its Igu apparent magnitude by applying K- 
corrections and colour terms based on a high-redshift single- 
burst model (zp = 5, T = 1 Gyr) which, as we have dis- 
cussed, reproduces the observed colours (Bower et al. 1992, 
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Figure 8. S(V — I) values for the HDF-N and HDF-S spheroidals com- 
pared with their clustered counterparts as a function 0/M450 absolute mag- 
nitude for three redshift ranges. Solid points represent cluster galaxies and 
open ones field spheroidals. 



Ellis et al. 1997). For the field spheroidals we computed M B 45o 
using B450 apparent magnitudes, given in version 2 of the 
HDF SExtractor catalogues, and a maximum likelihood 
method describe more fully in § 1.2, Anticipating this discus- 
sion, briefly the method allows us to obtain the most-likely 
K-correction for each individual pixel within the galaxy. This 
complication is necessary given the hybrid nature of the 



Figure 9. The observed 5(V — I) predicted for the blue-nucleated model 
spheroidal (solid line) compared with observed values in the HDFs. The 
colour gradient along the evolutionary track maps the observed evolution 
in y^06 — hli according to the right hand scale bar. Error bars in the field 
spheroidals show the lu deviation obtained from 100 bootstrap resamplings 
of the data pixels within each galaxy and field symbols follow the designa- 
tions of Figure 6. 



spectral energy distribution in the internally inhomogeneous 
cases. 

Clearly, over similar ranges of absolute magnitudes, 
field spheroidals show markedly larger dispersions in their 
internal colours than do cluster spheroidals. Although clus- 
ter spheroidals have more homogeneous internal colours, 
there is a trend for 5(V — I) to grow as a function of redshift. 
The trends calculated in Figure 2 shows that this is not an ar- 
tifact relating to a greater sensitivity to very small amounts of 
star formation. Unfortunately the limited depth of the cluster 
data does not allow us to put more meaningful constraints on 
the evolution of the cluster population. 



4 SPECTRAL SYNTHESIS MODELLING OF 
INTERNAL COLOUR VARIATIONS 

Simple comparisons between the internal colour distribu- 
tions of field and cluster samples allow us to establish, in a 
model-independent way, the proportion of active systems as 
a function of redshift and environment. Given that detailed 
spectroscopic information is not yet available for the sample, 
in order to learn more about the relevant timescales for the 
presumed star formation activity we must now turn to more 
detailed synthesis models. In this section we will investigate 
how simple models can be used to place broad constraints on 
the buildup of mass implied by the internal colour variations 
reported in this paper. 

4.1 A Simple Evolutionary Model For Field Spheroidals 

The blue light in field spheroidals reported in the present 
paper arises, in the majority of cases, in the central regions. 
Possible origins include residual star-formation from a recent 
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major merger (Kauffmann 1996), low-level star-formation 
that is a relic from a high-redshift collapse (Jimenez et al. 
1999), cooling flows from hot gas surrounding the central re- 
gion of the ellipticals, or a bursts of central star-formation 
possibly related to the inflow of gas arising from a recent mi- 
nor merger. In this section we aim to place constraints on the 
star-formation timescales and masses associated with the ac- 
tivity we claim to observe predominantly in the field popu- 
lation. We do this within the context of a simple model that 
seeks to reproduce the observed spatially resolved colours of 
the subset of active ellipticals with prominent blue nuclei. 

We assume that an intermediate redshift field spheroidal 
can be described with a two component model - an old un- 
derlying stellar component which formed at high redshift 
and a secondary younger component (responsible for the 
central blue light). The integrated colours of both compo- 
nents are estimated using the stellar population synthesis li- 
brary of Bruzual & Chariot (1996, BC96). We modelled the 
old component using a 1 Gyr exponentially declining burst of 
star formation commencing at a redshift of formation z F = 5. 
We adopt a flat Universe with O = 0.3, O a = 0.7. The sec- 
ondary population is modelled via an additional instanta- 
neous burst of star-formation presumed to occur close to the 
epoch of observati on. T he present modelling is similar to the 
one described in however here we include the evolu- 

tion as function of the galaxy age. 

In order to mimic the resolved properties of a model 
galaxy, we create an image of the two galaxy components 
using a modified version of IRAF MKOBJECTS. For the old 
component we use a de Vaucouleurs profile with a half-light 
radius of r e = 2.5 kpc. For the central blue component we 
also use a de Vaucouleurs profile, but with an effective ra- 
dius 0. lr e . We assume 15% of the galaxy mass is in the young 
component, and that the total mass is M = 1 x 10 12 Mq. 
Finally, between < z < z burst we combine the two compo- 
nents by adding up the images using the colours predicted 
by our spectral synthesis models. For the resulting galaxy we 
calculate 5(V — I) and repeat the procedure four times corre- 
sponding to 15% (by mass) bursts at redshifts z burst = 1.0, 0.8, 
0.6 and 0.4. 

Figure |9| shows the result of this exercise, superposed 
on the observed 5(V — I) for the HDF spheroidals. Error 
bars correspond to la estimates obtained by bootstrapping 
the data as described in Efron & Tibshirani (1986). We con- 
clude that our simple model can successfully reproduce the 
range of observed 5(V — I) as well as the observed increase 
5(V — I) with redshift. The essential point is that the ob- 
served internal colour scatter is not necessarily a relic of the 
initial elliptical formation event and can quite easily repre- 
sent a perturbation introduced by a recent episode of star- 
formation (forming only a modest proportion of the total 
stellar mass) superposed on a pre-existing old population. 
If this is the case then excursions from the homogeneous 
colours of the underlying old stellar population are proba- 
bly brief, lasting about 1 Gyr after which the galaxy returns 
to its pre-burst value. Ultimately high resolution spectra of 
large sample of spheroidals using Balmer lines will enable 
us to disentangle the internal clock of star formation for field 
spheroidals. However from this simple study it seems rea- 
sonable to conclude that as many as half the field spheroidal 
population share this activity by z ~ 1 . 



4.2 The Star Formation History of Spheroidals 

If our interpretation of Figure^] is correct, the question arises 
as to the extent of this recent activity in terms of the evolution 
of global star formation and its role for the mass assembly of 
field spheroidals as predicted in popular hierarchical models. 
Although the HDF samples are very modest ones on which 
to base such fundamental studies, we explore these wider 
issues using the analytical techniques described in Abraham 
et al. (1999). 

In the HDF-N, Abraham et al. were able to reconstruct 
a simple cosmic star-formation history by interpreting re- 
solved multi-colour data in the context of evolutionary syn- 
thesis models. For each HDF galaxy, an integrated star for- 
mation rate was determined and volume-averaged values 
(Psfr ) for the entire population was derived using a sim- 
ple V/V max methodology. Here we will attempt to estimate 
the contribution to Psfr arising from our HDF spheroidals 
and compare its evolutionary decline with time with that 
observed in other morphological populations. The crucial 
advantage of an analysis utilising the resolved colours is 
that it enables a determination of the relative ages and star- 
formations of the components within the galaxies. 

4.2.1 The Method 

Using the precepts of Bruzual & Chariot (1991), the total flux 
from a stellar system, F A observed at an age t can be de- 
scribed as the convolution of the spectrum of an evolving in- 
stantaneous burst, f\{i), with an star-formation history given 
by¥(t). 

Fa(to) = f'"y(t -t)f A {t)dt (3) 
Jo 

In our analysis, we use the stellar population synthesis li- 
braries from GISSEL96 (Bruzual & Chariot 1996) to estimate 
fx, assuming a Scalo IMF with lower and upper limits of 0.1 
and 125 Mq respectively and solar metallicity. For the SFR 
function W(t), we assume that the colours of a given stel- 
lar population can be modeled by exponential star-formation 
histories with a characteristic time-scale r of the form, 

V(t) = W e-'/\ (4) 

which enables us to describe a wide range of star-formations, 
given that t — > oo approximates a constant star formation, 
while t — > approximates an instantaneous burst. 

As the (J 30 o-band flux is weak for our sample, adding 
this to the analysis would enlarge the uncertainties in the ob- 
served colours of spheroidals rather than further constrain 
our estimates. Therefore, we concentrate on the B 450 , V 606 , 
and I$u bands only. We calculate the B 450 — V 606 vs Veoe — hu 
evolutionary tracks as observed at the galaxy redshift for 
a set of representatives formation histories (t = 0.1, 0.2, 
0.4, 1, 2, 4, and 9 Gyr) using the stellar population libraries 
from GISSEL96. For each spheroidal, we select all the contigu- 
ous pixels within a surface brightness threshold of B 450 
26 mag/arcsec 2 , and from the selected pixels construct the 
corresponding B i50 — V 606 vs V 606 — I 814 colour-colour dia- 
gram for each of the galaxies. We define this limit in order 
to maximize the signal/ noise associated per pixel in both the 
B450 — ^606 and V 60 6 — hn colours, selecting in the lowest sig- 
nal band in order to ensure an uniform signal criteria in the 
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Figure 10. A demonstration of our maximum likelihood methodology showing the observed and reconstructed Vgo6 — hu ana * B450 — Vgo6 colours. The 
reconstructed colours and SFR (lower right) maps were based on the the best fit models to the pixel colour-colour diagram (upper right). The colour assigned 
to the tracks drawn in the colour-colour diagram indicates age in Gyr according to the right hand colour bar (blue=young, red=old). As expected, the blue core 
is a region of current star formation, whereas the outer region shows virtually no activity. 



pixel selection. We estimate the optimum model and age for 
each pixel using a maximum likelihood estimator C defined 
as: 



C(t,r) 



2 

n 



ZnSQ 



exp 



2.5C, 2 



(5) 



where we take the product over both colours B i50 — Vgo6 
and V 606 — J 814 , C, and 5C, represent the data colour and 
colour error associated with the pixel respectively, and c, is 
the colour of the modelled track at a given age and e-folding 
time-scale T. For each pixel we obtain the most likely age, 
SFR and t value. Repeating over all acceptable pixels, we ob- 
tain the integrated value at the observed redshift. This pro- 
cedure is repeated for all spheroidals in the HDFs. In order 
to illustrate this methodology, in Figure [lo] we show a typi- 
cal pixel-by-pixel colour-colour diagram and illustrate how 
the predicted evolutionary tracks at the redshift of observa- 
tion can reconstruct the V 606 — 7 8 i 4 and B 450 — V 606 observed 
colours and construct a SFR map for each galaxy pixel. 

Given that our methodology depends explicitly on the 
redshift of the galaxy, we need to take into account the uncer- 
tainties associated with the photometric estimates when cal- 
culating SFR for galaxies with photometric redshifts. Up to 
z ~ 1, photometric redshifts have proved to be in reasonable 
good agreement with spectroscopic redshifts with an rms er- 
ror of <5z ~ 0.1 for red populations (Hogg et al. 1998, Gwyn 



et al. 2000, de Soto et al. 1997). We performed Monte Carlo 
simulations by refitting each galaxy using the models for the 
appropriate photometric redshift plus a random component 
of 5z = 0.1. Repeating this exercise for each galaxy with pho- 
tometric estimates allows us to obtain an uncertainty on the 
total SFR estimated for those systems for which only photo- 
metric redshifts are available. 

Clearly the predictions of spectral synthesis models are 
sensitive to both metallicity and dust content, neither of 
which is well constrained in our data. The presence of dust 
would have the effect of underestimating the computed SFR, 
whereas metallicity changes will produce systematic shifts 
between the predicted age and colours of each galaxy. Whilst 
recognising the limitations involved, our primary goal at this 
stage is to illustrate the potential of our resolved colour tech- 
nique to determine redshift-dependent trends rather than to 
derive precise absolute values. Throughout the SFR analy- 
sis we will use solar metallicity in our calculations, as from 
Figure ^ we can see that this reproduces fairly well the ob- 
served colour of low 5(V — I) spheroidals. With resolved in- 
frared data such as may be available from NICMOS, it may 
ultimately be possible to limit the uncertainties arising from 
dust. 
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Figure 11. TTie comoving SFR density computed using the 63 spheroidal 
(filled circles) and 51 morphologically-irregular (filled stars) galaxies in 
both HDFs. Corrections for incompleteness were made assuming SFR oc L, 
local luminosity functions and evolutionary corrections described in the 
text. Hollow symbols represent uncorrected values. Error bars take into ac- 
count systematic uncertainties implied in the calculation of the SFR in in- 
dividual galaxies and bootstrap errors which account for the small sample 
size. 



4.2.2 Evolution in the Star-Formation density of spheroidals 

We now estimate the SFR density (psfr) of spheroidals as 
a function of redshift and choose to compare these esti- 
mates with similar ones derived for a parallel sample of 
morphologically-irregular galaxies in both HDFs. In the case 
of the latter, we used the morphological classification of 
Abraham et al. and adopted the same Igu limit for both 
HDFs. 

We calculate the SFR density at a given redshift interval 
Z\ < z < z 2 using the expression, 



Psfr(z) = £ 



\)Z\<Z{<Zi 



SFR, 

VmaxA 



(6) 



where the sum is over all the galaxies i for which zi < z, < z 2 
and V max is the volume at z max , the largest redshift at which 
the galaxy i remains brighter than our limiting magnitude 
(Schmidt 1968). Clearly as we move towards higher redshifts 
our sample becomes increasingly biased towards intrinsi- 
cally more luminous galaxies and we will consequently ne- 
glect a contribution to Psfr from sub-luminous objects. To 
correct for this we need to assume: i) the dependence of the 
SFR on the luminosity L, ii) the shape of the spheroidal lu- 
minosity function and Hi) the luminosity evolution assumed 
for the spheroidal population. Whilst none of these is un- 
fortunately known precisely, we can make progress by as- 
suming that SFR oc L (i.e. the most luminous galaxies form 
the most stars) and adopting a suitable Schechter luminosity 
function (LF). Under these assumptions, the corrected SFR 
density (psfr) would become: 



Psfr(z) = C(z)psfr(z) 



with, 

C(z) = 



Jo 00 SFR(L)qb(L)dL _T(a + 2) 



/" SFR(L)cp(L)dL 



(7) 



(8) 



10 0.4(M*-M) (a+ 2) e -10°«« -M) dM 

-00 



(9) 



where M* may be permitted to evolve with redshift (e.g. 
according to an evolutionary correction) and M~ and L m ;„ 
are the faintest absolute magnitudes and luminosities respec- 
tively within the redshift interval in question. 

As an illustration of the technique, we adopted 
Schechter LFs from Marzke et al. (1998) with M* b . = -20.87 
and a = — 1 for the HDF spheroidals and M ( * = —21.28 and 
oc = —1.81 for the irregulars. For the spheroidals we as- 
sumed M* evolves according to that appropriate for a pas- 
sively evolving system with (t = 1 Gyr). In the case of the 
irregulars, we assumed a constant SFR (t = 5 Gyr). 

Finally is necessary to define a number density correc- 
tion factor, /C(z), very similar to the one described in equa- 
tion to estimate the redshift completeness of our sample, 



K(z) = 
and 

B = 0.4 



r(ce + l) 
B 



ln(10) / 

J — C 



10' 



0.4(M*-M)(a+1) 10"- 4 (M*-M) 



dM 



(10) 



(11) 



where M~ is the faintest absolute magnitude included in the 
sample at a given redshift determined by our flux limit of 
Igi4 = 24 mag and the distance modulus relation. 

Although /C(z) and C(z)clearly depend on the poorly- 
constrained LF shape, for our adopted values we find com- 
pletenesses (in terms of the contribution of detectable galax- 
ies to the integrated p SFR ) of ~ 90% and ~ 80% for 63 HDF 
spheroidal and 51 HDF morphologically-irregular galaxies 
respectively. 

Figure [ll] shows the results of our calculations and the 
uncertainties introduced by our corrections for incomplete- 
ness. As expected these corrections become more impor- 
tant at higher redshifts and for the morphologically-irregular 
population whose LF is presumed to be steep at all times 
(a < —1). However in most of the cases the corrections are 
comparable to the size of the statistical uncertainties. Ran- 
dom errors bars were calculated by adding in quadrature the 
systematic errors in the SFR calculation for individual galax- 
ies, including the Monte Carlo estimates for the photometric 
redshift and Poisson errors estimated by bootstrapping 100 
times the sample when computing the SFR density. In Ta- 
ble | we list the computed values for p SF R as well as their 
uncertainties. 

Figure ^ enables us to provide the first rough con- 
straints on the contribution of the increased colour variation 
in the HDF spheroidals to the evolution of the 5(V — I). A 
rise with redshift is discernible but the uncertainties are con- 
siderable. Interestingly, the evolution in psfr mirrors, but at 
a lower level, that established more convincingly for field ir- 
regulars from the CFRS/LDSS sample of Brinchmann et al. 
(1998) (but estimated here for the HDF from a completely in- 
dependent resolved colour method). Brinchmann et al. com- 
puted p SFR between 0.3 < z < 0.9 using the equivalent width 
of [Oil]. Although their p SF R values are systematically lower 
for both the spheroidal and peculiar /irregular population, 
similar evolutionary trends were observed in the overlap- 
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ping redshift range 0.4 < z < 1.0. If it could be established 
that the decline in the spheroidal p SFR was driven by that as- 
sociated with the demise of the more active irregular popu- 
lation, this would provide further support that the blue cores 
observed in the HDF spheroidals arise from recent mergers 
of the irregular population (Le Fevre et al. 1999, Brinchmann 
1999, Brinchmann & Ellis 2000). 



activity, it may be possible to strengthen the connection 
between the declining luminosity density contributed by 
morphologically-irregular galaxies and the implied contin- 
ued assembly of field ellipticals. The goal remains an impor- 
tant one in understanding the origin of the Hubble sequence 
and in defining the role of the environment (Brinchmann & 
Ellis 2000). 



5 CONCLUSIONS 

We have exploited HST's unique capabilities to resolve dis- 
tant galaxies by examining the internal colour distribution of 
a sample of 79 field spheroidals in both Hubble Deep Fields. 
From a model-independent analysis of our data, we find: 

• A significant fraction (>30%) of field spheroidals show 
marked variations in their internal colours. In most of the 
cases where colour variations are seen, this manifests itself 
in terms of centrally-located blue cores. 

• Using a statistic, 5(V — I), to characterise the internal 
homogeneity of a galaxy we have extended our analysis to 
permit a comparison of the internal variations found in the 
HDF sample with that for a sample of spheroidal galaxies in 
five rich clusters spanning a similar redshift range. Marked 
differences are found in the sense that galaxies of the same lu- 
minosity and redshift are consistently more uniform in rich 
clusters than in the field. Assuming the blue colours located 
in the HDF field ellipticals arise from recent star formation, 
our comparison provides strong evidence for more extended 
star formation histories in field ellipticals (as expected in hi- 
erarchical models). 

• We compare our observational diagnostics with those 
used by other workers and find general agreement. The most 
straightforward interpretation of our data is that a fraction 
of field spheroidals at intermediate redshift has undergone 
recent star formation, perhaps as a result of recent merging 
or an inflow of material. 

Using spectral synthesis models and a maximum like- 
lihood analysis developed by Abraham et al. (1999) we at- 
tempt to reproduced the observed colour variations for the 
field spheroidals and to compute their current SFR. From this 
model-dependent analysis we find that: 

• Assuming that the observed blue nuclei observed in 
some field spheroidals arises from a recent burst of star- 
formation superposed on a pre-existing old stellar popula- 
tion, we can readily reproduce the observed distribution of 
6(V — 7) values. The star formation episodes implied are 
short, lasting less that 1 Gyr, and typically involve 15% of 
the galactic mass. 

• If the HDF fields are representative, the fraction of field 
ellipticals undergoing activity of this kind could be as high as 
~50%atz~l. 

• The contribution to the global SFR density from field 
spheroidals, whilst uncertain to determine accurately be- 
cause of the necessary corrections for incompleteness and 
possible luminosity evolution, appears to rises modestly as a 
function of redshift in a manner which may connect with the 
demise in activity associated with morphologically-irregular 
systems. 

With a larger sample, augmented with diagnostic spec- 
troscopy capable of constraining the timescale of recent 
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Table 2. SFR density estimates 
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1.05 -2.106-jll -2.0701^1 -1.442:} : °J° -0.883I? : « 2 



Single Burst (t = 1 Gyr), z F = 5, Marzke et al. E/SO LF (as = -1, M* = -20.87) 
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